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Abstract 



different methods. In the spirit of a many-body approach we make an expan- 



The temperature evolution of the quark condensate is studied using three 

O 
l> . 

sion in the scalar density up to second order. Our result is consistent with 
\ chiral perturbation theory to two-loop order. 



1. INTRODUCTION 



Chiral symmetry, which is an approximate symmetry of QCD in the light-flavor sector, 

X 

is spontaneously broken in the QCD vacuum. A topic of interest for nuclear physics is the 
restoration of chiral symmetry, at high temperature and/or large baryonic density. The 
approach to restoration shows up as a decrease of the magnitude of the quark condensate, 
which is the order parameter of the spontaneous breaking. Within a hadronic picture the 
physical origin of this restoration is very simple in nature. Each hadron contributes a 
certain amount to the decrease, governed by its sigma commutator E^. For an assembly 
of independent hadrons, their effects add. This leads to the following expression for the 
evolution with temperature of the condensate 

h -ITT " t 7T 



where the sum extends over all hadron species present in the medium and p s h is the cor- 
responding scalar density. In a baryon free regime (zero baryonic chemical potential) the 
scalar density is produced by thermally excited hadrons. At low temperature eq. (|1]) is the 
leading contribution. Deviations arise from the interaction between the hadrons and the 
major problem is to evaluate their effect in a way which is consistent with chiral symmetry. 
In this work we shall investigate the temperature evolution of several quantities linked to 
chiral symmetry, namely the quark condensate, the pion mass and the weak pion decay con- 
stant. The results that we obtain in the chiral Lagrangian approach have already been given 
in chiral perturbation theory [|l],|2|. Our approach provides alternative methods suitable for 
incorporating other hadron species, especially vector mesons, in a consistent way. 

For temperatures less than the pion mass the presence of mesons other than pions can 
be ignored, since their thermal abundance is small. It is also sufficient to expand to second 
order in the density. 

To respect chiral symmetry we start from the standard nonlinear chiral Lagrangian, 
including the explicit chiral symmetry breaking piece, namely 

f 2 1 
£ = J -± Trd^Ud^ + -flml Ti(U + U ] ). (2) 

The 5*?7(2)-matrix U is usually expressed in terms of the canonical pion field $ as 

U = exp (ir.^F(X)) (3) 

where F(X) is an arbitrary odd polynomial of X = ($ 2 // 2 ) 1/ ' 2 

F(X) = X + aX 3 + (3X 5 + (4) 

Different choices of the function F(X) have been introduced in the literature, corresponding 
to different values of the parameters a and (3. For instance, the simple choice F(X) = X 
is commonly used in chiral perturbation expansions. The choice of Weinberg is F(X) = 
arcsin[X/(l + X 2 /4)] §. On the other hand F{X) = arcsinX gives the PCAC field, 
satisfying d^A? = —m n f%&i. Physical observables should not depend on the particular 
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function F(X). In practice, however, a particular choice is made and some approximations 
are needed. It is thus important to check the independence of the results on F(X). 

In the case of a cold baryonic medium the independence on the representation has been 
established recently by Delorme et al. j|] and by Kirchbach and Wirzba ||. It is one of the 
purposes of the present paper to explicitly check this property for the quark condensate in a 
hot baryon-free environment. A baryon-free zone may possibly be realised experimentally at 
RHIC and the LHC and is certainly present in the early universe. To derive the temperature 
dependence of the quark condensate we will use three methods. The first (sect.l) relies on 
the grand potential method introduced in ref. ||. In the second method (sect. 2) we evaluate 
the thermal expectation value of the chiral symmetry breaking piece of the Lagrangian. 
Both methods demonstrate the expected independence of the result on the particular choice 
of F(X). Once this independence is proven we will introduce a third method (sect. 3) which 
requires the validity of the PCAC relation, i.e. corresponds to a particular choice of F(X). 
In this case the condensate can be derived as the scattering amplitude of soft pions in the 
system. The condensate is then expressed in terms of the self energy for soft pions. In this 
derivation we will prove the validity of the Gell-Mann-Oakes-Renner relation (GOR) in a 
hot medium, extending the previous derivation of Chanfray et al. |7|] which was restricted 
to zero temperature. 



2. THE GRAND POTENTIAL METHOD 

To obtain an expression for the chiral condensate ratio, (qq)T/(qq)oi i n terms of pionic 
degrees of freedom we can use the fact that, according to the Feynman-Hellmann theorem, 
(qq)r is related to the derivative of the free-energy density (the grand potential) £l(T) with 
respect to the bare quark mass m as 

(qq) T = dQ(T)/dm (5) 

(we ignore the fact that up and down quark masses not the same). Denoting the difference in 
free energy density as f2(T) = &(T) — Cl(0) and using the Gell-Mann Oakes Renner relation 
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(GOR) in the vacuum 



m lfl = -2m(gg) = ~ m ~^- ( 6 ) 



one immediately derives that 

(qq)o f 2 dml 

This expression is identical to that of ref. if one ignores small corrections to f n and m n 
of order m. For a free pion gas the grand potential is 

n(T) = — j Q dk k 2 In (l-e-^ T ); u k = ^ + k*. (8) 

Taking the derivative with respect to the pion mass and using ([7|) leads to 

(qq)t = 1 _ m n /2 
(qq)o flm\ 

with 



PliT) (9) 



pliT) ~ 3 / (2^)3 Zl expK/T) - 1 • (10) 
Keeping in mind that the pion sigma commutator is = (with nonrelativistic nor- 

malization) and that p s n is the scalar density of pions, we recover the independent-particle 
expression of eq. (|J) for a gas of free pions. 

In the evaluation of the free energy, Q(T), for interacting pions we start from the well- 
known expression for interacting boson systems at finite temperature fl9",IQfl. In terms of the 
pion Green's function 

D(lu, k, T) = (cu 2 -ml-k 2 - U{u, k, T))-\ (11) 

the grand potential is given by 

fl(T) = tf(T) -Uj^J^ m Imjlnl-D- 1 ^, k, T)] + D(u, k, T)U(u, k, T)} (12) 

where f(uj) = (exp(u/T) — 1) 1 is the Bose factor and fl'{T) is the sum of all contributions 
from 'skeleton diagrams' representing the perturbation expansion of Q. These are evaluated 
by using full single-particle Green's functions rather than bare ones. Furthermore 



sn/su = o 



(13) 



and 



sn'/5D = n 



(14) 



(see refs. [PHTTH ) . Since the physical pion mass enters the propagator, D, the self energy 
n(u;,k, T) is the difference 



IL(u, k, T) = n(w, k, T) - n(w, k, 0), 



(15) 



whereby the infinite contributions to the mass operator H(u>, k, 0) due to vacuum fluctuations 
are removed. Consequently also fl{T) is finite and no regularization of the loop integrals is 
needed. 

The key quantity in eq. (0) is the pion self energy EL To calculate it we use the nonlinear 
ex-Model Lagrangian defined in eq. (2). To lowest order, the interaction part takes the form 
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ml [a- — ] ($.$) 2 + ( a - - ) M <^$.<9 M $ + ( 2a + - ) $.<9 M $ 
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(16) 



The pion self energy EI contains single-pion and three-pion contributions (see Fig. la and 
lb) 



IL(u, k, T) = n(w, k, T) + n(w, k, T) d7r . 



(17) 



In the present approach and for the temperature range of interest (T < m^), only the 
single-pion contribution is needed. This is equivalent to the Hartree approximation and 
leads to : 



U(co,k,T) = Xj(T) 



10\ 9 
]Qa - —j ml 



lOd - - ) (ml + LU 2 -k 2 " 



(18) 



where A = 1/Qf 2 and j(T) is the thermal loop integral given by 

dk 1 1 



m = 3 / 



(27r) 3 2u k exp(u k /T) 
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(19) 



Here ui^ = (to 2 -l-k 2 ) 1 / 2 is the in-medium pion dispersion relation with an effective mass m n . 
Its value is obtained by solving the following equation : 

ml = m 2 +U(m 7r , k = 0,T) (20) 

which yields 

^ _ 1 + (80a - 10/3)Aj m2 
n 1 + (80a - 16/3)Aj n V ' 

Since m 2 is proportional to the bare mass it vanishes in the chiral limit m n — > 0, at all 

temperatures as required by the Goldstone theorem. Notice that, for 'on-shell' pions, the a- 

dependence of the self energy disappears, as it should, and U os = 2Xj. The pion propagator 

D(u,h, T) can now be written as 

D(^,k,T) = 7 (^ 2 -k 2 -m 2 )- 1 , (22) 

In addition to the effective mass it also contains a temperature dependent residue || 



7(T) 



40a - jj) Xj 



-i 



(23) 



which arises from the energy dependence of the 7T7r-interaction given in (|1G|). 
Finally we need to evaluate the interaction part Q' which takes the form 
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x ( 80a - z —\mt- (40a 



l—QO 

y-k 2 )) + (^-p 2 )) 



(24) 



3 7 * V 3, 

Once £l(T) is determined one can find the chiral condensate ratio from eq. (^). When taking 
the derivative one should note that Q depends on the free pion mass m % , both explicitly 
and through the self energy IT. Due to the stationarity conditions (p~3| ) and ([14]) the in- 
dependence of II does not influence <9f2/<9m 2 and we obtain 

= d( T) - ^(80a - 10/3)d 2 (T). (25) 

where d(T) = ~f(T)j(T) which is identical to the result of ref. H for a = —1/12 (the 
'Weinberg' choice for F(X) 0). Hence, with eq. (^), 



= 1 - QXd(T) - (240a - 10)(Ad(T)) 2 (26) 

At first sight, this expression has an unwanted a- dependence. This model dependence 
disappears, however, if we expand the residue given in eq. (|23"D to first order in the coupling 
constant A. This gives the evolution of the condensate to second order in A, independent of 

a 

M^ = 1-6Aj(T)-6(Aj(T)) 2 (27) 
\QQ)o 

where the first-order term 6Aj is just the contribution from a free pion gas (see eq. @), 
which dominates. The expression (|27|) represents, in the spirit of a many-body approach, 
an expansion in the density. We can check its consistency with chiral perturbation theory. 
In the chiral limit (m n — > 0) where j can be calculated analytically (j = T 2 /8), we recover 
the chiral perturbation result to order T 4 

fah _ ! T> T 4 

(gq)o W 384/^ 1 ) 

of Gasser and Leutwyler |[|. 

3. EXPECTATION VALUE OF THE CHIRAL SYMMETRY BREAKING 

HAMILTONIAN 

The condensate ratio can also be obtained directly from the in-medium expectation value 
of the symmetry breaking Hamiltonian of the nonlinear sigma model (eq. (2)) 

H xSB = -\mlfMU + tf) = -ml fl cos F{X). (29) 

Using the GOR in vacuum (eq. ((!])) we obtain 

M^ = (cos(F(X))) T (30) 

where {..)t on the rhs denotes a thermal expectation value. In order to be consistent with 
the second-order expansion in expression (27) we have to evaluate the thermal expectation 



value to fourth order in the pion field : 



(Wh , 1,4> 2 , / 1\ 4> 2 * 2 

In order to evaluate the expectation values explicitly, the pion field is expanded in terms 
of creation and annihilation operators Bj k and Bi ^ for the quasi pions, specified by the 
propagator given in (f2^) 

0=7'/ ~ ^ , 1/2 (^k(*) + Sj_ k W) (32) 

The presence of the 7 1//2 factor ensures that the form ( p2|) of the pion propagator Z) is 
recovered once the time dependence B^(t) = B^ exp(— iu^t) is introduced. We have verified 
that the field operator (|32]) obeys canonical commutation relations. The expectation value 
of ($ 2 )t is calculated by normal ordering with respect to physical vacuum and yields 

(^^/p^-E^Wr^^ (33) 

Similarly ($ 2 $ 2 )r is obtained by using Wick's theorem and gives 

($ 2 $ 2 ) T = jj ($ 2 )| = ^ 7 2 j 2 (34) 



Expanding the residue 7, defined in (j23"|), for both ( |33| ) and ( p4|) to first order in Xj and 
inserting in (|3"1"D, we recover the expression Q27| ) of the previous section. In the next section 
we will introduce a third method which is based on the use of the PCAC pion field. 

4. SOFT-PION METHOD 

In this section, the values of the parameters a and (3 in the expansion of F(X) (eq. (|4])) 
are chosen such that the PCAC relation 

d ll At(x) = -Uml<$> i {x) (35) 

is fulfilled. This corresponds to the values a = 1/6 and (3 = 3/40. In the previous sections, 
however, it has been shown that the condensate ratio, to second order in the density, is 
independent on these parameters, which legitimates this third approach. The aim here is to 
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obtain the condensate ratio from its relation to the self energy for soft pions. In deriving such 
a relation we will show the validity of the in-medium version of the GOR which relates the 
effective pion mass, the effective pion decay constant and the in-medium quark condensate. 
This has been previously shown to hold in a dense but cold hadronic medium by Chanfray 
et al. 0. 

We start from the QCD operator identity 



Q, 5 (0),^f(0) 



-i2mqq(0). (36) 



Taking the thermal expectation value of this relation and using closure on the lhs one has 
~Z 



dv Y, e —^- ((n|A (r,0)|m}(m|^(0)|n) - (n|^(0)|m}(m|A (r, 0)|n}) (37) 



n,m 



= -2im(qq) T (38) 

where Z is the partition function. After spatial integration, translational invariance requires 
that p n = p m . 

The above equation is equivalent to the following energy-weighted sum rule 

2m(qq) T = - £ -— 2(E m - E n ) |(m|4°(0)|n)| 2 , (39) 

n,m 

It is now possible to derive a generalized GOR relation at finite temperature. For this 
purpose we assume that the sum rule (|39| ) is saturated by exciting a single quasi pion, n, 
(of mass m n ) from the heat bath, i.e. 

g — /3E n 

2m(qq) T = - ^ —— 2 rh n | < n + 7f| A$(0) \n > \ 2 (40) 
It is natural to introduce a thermally averaged pion decay constant, f n , through 

E^l<- + -I4°(0)|n>| 2 = iZgL (4i) 

The generalized GOR relation follows immediately 

2m(qq) T = ~]lm\. (42) 



Notice that Lorentz invariance is broken by the very concept of temperature. Therefore the 
renormalization of f n from the GOR relation ( JI^ ) concerns only the time component of the 
axial current. 



Coming back to the energy- weighted sum rule (p9|), we introduce the PCAC interpolating 
pion field through = —d fJi Af(x)/f n m'^. Taking the matrix element of this PCAC 

identity between states \m) and |n) yields 

t(E n - E m ) (n\J$(0)\m) = -Uml (n|^(0)|m) (43) 

and leads to an alternative form of the sum rule 

2m (qq) T = -2fim* ]T — — ' ^ " . (44) 

This relation links the condensate to the axial polarizability of the thermal bath. On the 
other hand it is a simple matter to show that the thermal pion propagator Z}(0,0, T), 
calculated at the soft pion point u = 0, k = 0, is 

£>(o,o, t) = -K + n^o.T))- 1 

dt f dT({-i)[$i(T,t), $i(O)]) T 0(*) 



_ 2 y^^^NMM (45) 

n.m ^ m -^m -^n 



Comparing eqs. ( 44j) and (ffq), we obtain 



<*>r ^ + n ( o,o r)y 



(?9>o V , 
This is the generalization of a result already obtained for nuclear matter by M. Ericson 



I4 , which incorporates the coherent rescattering of the soft pion in an optical potential 
coinciding with the soft-pion self energy 11(0, 0, T). 

In order to evaluate the condensate, we need the expression of the self energy to two-loop 
order, including in particular the contribution from the diagram of Fig. lc. These were not 
necessary in our previous approach ([18]). The two- loop contribution requires to expand the 
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Lagrangian to sixth order in the pion field. Both parameters a and (5 enter in this case. The 
result for the self-energy is 

- Q(X 2 ) r + ^(X 2 )|)(^-k 2 ) (47) 

where (X 2 )t = 12Aj represents the thermal loop integral of quasi pions. In expression ( f47D 
we take for consistency the effective mass to first order, i.e. m 2 = m 2 (1 + 2Xj). We obtain 
the self energy for soft pions from (f47|) and when inserting its value in ([|6|) we recover, 



to second-order in Xj, our previous result of eq. ([H]). Notice that the expression (|47|) for 
n(o;, k, T) is valid at the level of an extended Hartree approximation. To this order there is 
an additional contribution corresponding to a pion fluctuating in three pions (see Fig. lb) 
which, however vanishes at the soft pion point. Therefore it does not affect the soft-pion 
self energy and hence the condensate. However it enters in the residue 7 to second order in 
the coupling constant A. Since we did not calculate this term, our residue and hence our 
expression (21) for the pion mass are only valid to first order. This is also the case for the 
pion decay constant From the GOR relation we obtain f n = f T (l — 4Aj), in agreement 
with ref. 0. 

In summary we have studied the evolution with temperature of the quark condensate in 
the framework of a many-body approach, considering thermal excitations of pions. The nir 
interaction is taken from the nonlinear sigma model. Our expansion parameter is the scalar 
density of the quasi pions and we work up to second order. We ave shown in two different 
approaches that, up to this order, the condensate is independent of the transformation of the 



pion field as expected. In line with ref. |T2j we have used a third method, based on PCAC, 
to derive the condensate which is related to the self energy of soft pions. In the chiral limit 
our results are consistent with chiral perturbation theory, to two-loop order. The next stage 
will be to introduce the rho meson exchange in the tttt interaction in a way which preserves 
the chiral properties of this interaction. 
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Figure Captions 



Fig. 1: The pion self energy including the tadpole (a) 3-ir intermediate states (b) and the 
two-7r loop contribution (c). 
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FIGURES 
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